I 



FLYBY MEASUREMENTS OF THE PROFILES OF 
WATER VAPOR, PARTICULATE MATTER, AND 
CARBON DIOXIDE IN THE MARTIAN ATMOSPHERE 


by Alfred E. Barrington, Anthony J. Caruso, 
and Gene G. Mannella 

Electronics Research Center <;? 

Cambridge, Mass. . 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION • WASHINGTON, D. C. AUGUST 1967 


TECH LIBRARY KAFB, NM 



TECH LIBRARY KAFB, NM 



FLYBY MEASUREMENTS OF THE PROFILES OF 
WATER VAPOR, PARTICULATE MATTER, AND 
CARBON DIOXIDE IN THE MARTIAN ATMOSPHERE 

By Alfred E. Barrington, Anthony J. Caruso, and Gene G. Mannella 

Electronics Research Center 
Cambridge, Mass. 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


For sale by the Clearinghouse for Federal Scientific and Technical Information 
Springfield, Virginia 22151 — CFSTI price $3.00 


I I I 




FLYBY MEASUREMENTS OF THE PROFILES OF 
WATER VAPOR, PARTICULATE MATTER, AND 
CARBON DIOXIDE IN THE MARTIAN ATMOSPHERE 


By Alfred E. Barrington, Anthony J. Caruso, and Gene G. Mannella 

Electronics Research Center 


SUMMARY 

A study of infrared absorption measurements in the Martian atmosphere for the 
determination of the profiles of water vapor and carbon dioxide is presented in quanti- 
tative detail. 

Using the best available data on the Martian atmosphere, it is shown that with 
state-of-the-art instrumentation the following can be measured: an upper limit of 
water vapor of 10 -3 g/cm2-column (at 6.3 /i), scattering by particulate matter 
(at 3. 5 fj.), and the carbon dioxide profile above an altitude of 50 km (at 4.3 /u). 

The choice of wavelengths and of experimental instrumentation is based on a realis- 
tic appreciation of Martian atmospheric characteristics and experience gained during 
the Mariner II Venus mission. 


INTRODUCTION 


The most exciting aspect of Martian exploration is the search for extraterrestrial 
life. The prerequisite for the formation and sustenance of living organisms is water. 

So far, the only acceptable evidence of water on Mars is the observation by Kaplan et 
al (ref. 1) who, with a terrestrial telescope, were able to record rotational absorption 
lines characteristic of water vapor near 0. 82/Lt. From their observations they derived 
a figure of approximately 10 “3 g/cm^-column of water vapor in the Martian atmosphere. 


It was assumed until recently that the Martian atmosphere consisted mainly of 
nitrogen. However, the Mariner IV observations of the refractivity and density of the 
Martian atmosphere have led to the conclusions that carbon dioxide is by far the most 
abundant constituent (ref. 2). In constructing a Martian model atmosphere based on 
Mariner IV measurements, Johnson (ref. 3) has indicated an inconsistence be- 
tween the probable temperature distribution through the Martian atmosphere and the 
above mentioned spectroscopic determination of water vapor on Mars. The question 
regarding the water vapor content of the Martian atmosphere is thus of paramount 
interest to planetary physics as well as to extraterrestrial biology. The earliest 



opportunity for a measurement in the vicinity of the planet will probably occur during a 
flyby. The present study was initiated with this purpose in mind. 

With the feasibility of a flyby experiment established, a detailed proposal was sub- 
mitted by the authors to the Office of Space Science and Applications in March 1966. * 

In essence, it was proposed to measure the absorption of solar infrared radiation by 
water vapor (6. 3 jz) and also to obtain a simultaneous measurement of the carbon diox- 
ide profile (4. 3 jx) in order to corroborate the interpretations of the Mariner IV oc- 
cultation experiment. A suggestion by Ben-Aryeh (ref. 4) to perform a similar meas- 
urement for water vapor at 6.3 /j. from a Martian orbiter has recently appeared in the 
literature (September 1966). He also suggested the need for an absorption measurement 
at 8.4/z, where the Martian atmosphere should be completely transparent, to determine 
the effects of scattering on absorption measurements in this portion of the infrared 
spectrum. However, his discussion of absorption measurements is couched in only 
very general terms, whereas the authors' earlier study contained considerable quanti- 
tative detail. This shows, for example, that, in the infrared, Rayleigh scattering is 
insignificant in the Martian atmosphere. Only particulate scattering, as perhaps from 
frozen carbon dioxide or dust particles might be observable. Furthermore, solar radi- 
ation at 8. 4 /lx is insufficient for meaningful measurements with state-of-the-art detec- 
tors. For adequate intensity, radiation in the atmospheric window near 3. 5 /n should be 
used instead. 

The quantitative evaluation of infrared absorption measurements in the vicinity of 
Mars as presented here is considered of sufficient interest to warrant distribution 
among the scientific community. 


THEORY 

Water vapor and carbon dioxide have strong infrared absorption bands. Therefore, 
the profiles of these constituents of the Martian atmosphere can be determined by ap- 
propriate measurements of absorption of solar radiation during periods immediately 
prior to or following occultation of the sun by the planet. In this case, absorption of 
solar radiation will occur along a path parallel with a tangent to the planetary disc and 
will be a maximum along the tangential path. 

The following calculation of the pertinent infrared absorption characteristics of the 
Martian atmosphere is based on the barometric equation and the best available figures 
for surface number density, water vapor content, Martian atmospheric temperature, 
and other well established Martian constants. 


♦Water Vapor and Carbon Dioxide Profiles of the Martian Atmosphere, NASA-ERC 
proposal, 14 March 1966. 
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Let n (h) = number density at altitude h 

n * number density at surface (altitude = 0 ) 
s 

H = scale height 

“16 

k = L 38x10 erg/°K 

T = average atmospheric temperature 

m = mass of molecular particle 

g Q = surface value of gravitational acceleration. 

Then (see ref. 5): 

h_ 

n(h) = ne H (1) 

s 

where 



( 2 ) 


2 

The number of molecules Njj/cm -column above the surface is obtained by integrating 
equation (1), i.e. : 



o 


The mass Mjj/cm^-column above the surface is given by 

M = n mH . (4) 

H s ' ' 

2 

The number of molecules N z and the mass M z /cm -column along a path z through the 
Martian atmosphere, parallel with a tangent to the planetary disc (radius R) and spaced 
from it by a distance h Q , is found from Figure 1 as follows: 
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Figure 1. 


z 2 = (R + h) 2 - (R + h Q ) 2 2 2R(h - h Q ) for h «R, 


or 

h 


= h 


o 



(5) 


The path extends from z = - ro toz= + C0 . 
Substituting equation (5) in equation (1): 


n(z) = n g e 


|-h 

o 

_ H 


2 1 
z 

2RH. 


and 


N = n 
z s 





( 6 ) 


(7) 


Let 

= t 2 

2RH 

Therefore, dz = V2RH dt 

Substituting equations (8) and (8a) in equation (7) : 

h 

o 

N = n V2RH e H 
z s 



( 8 ) 

(8a) 


(9) 
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Maximum absorption occurs when the path is along the tangent, i. e. , when 


h = 0 . 
o 


Equation (9 ) then reduces to 


N = n \2 ttRH = 2. 5n Vrh" 

z, max s s 

M = 2. 5 mn Vrh 

z, max s 


(9a) 

( 10 ) 


In order to utilize published values of infrared absorption of water vapor and car- 
bon dioxide, we can use the following information (see ref. 2): 

"3 2 

For water vapor: M = 10 g/cm - column 

ri 

17 3 

For carbon dioxide: n ^ 2 x 10 mol/cm 

s 

Surface pressure: P g ~ 5 mb 


Scale height: 


H ~ 10 km 
C °2 


Martian radius: 


R = 3360 km 


Absorption by Water Vapor 

From equation (4), for water vapor: 

.-3 


mn = 
s 


10 


10 x 10 


i.-9 . 3 

= 10 g/cm . 


Substituting this value in equation (10): 

M = 2. 5 x 10" 9 \3360 x 10 x 10 

z, max 


10 ) 


1/2 


= 0.046 g/cm -column (11) 



This latter quantity is also the amount of precipitable water, W (pr-cm). The 
transmittance T^ at 6. 3 ji of precipitable water vapor can be calculated according to 
the following expression (see ref. 6): 

V -1. 97 W/W q / 

“ exp \ r — ' i/2° 1/21 1/2 t (12) 

X I [l + 6. 57 (P q /P) (T/T q ) 1/2 (W/W o ) 1/2 J 1/2 J 

where 

W = water vapor concentration for 50% transmittance at 6. 3/i 
0 (0.002 pr-cm) 

at 

P = 740 mm Hg 
o 

and at 

T = 273° K . 
o 

According to Johnson (ref. 3), there is little, if any, water vapor above 20 km. We 
assume an average pressure P of 1 mm Hg and an average temperature T of 180° K 
(ref. 2) over the altitude range of interest. 

Using these values of P and T and the value of W from equation (11), we obtain: 

W/W = 23 
o 

P /P = 740 
o 

T/T =0.66 . 
o 


Upon substituting these quantities in equation (12), we obtain: 

-45.5 


T^ = exp 


= exp 


_(1 + 6 . i 


57 x 740 x 0.814x4.8) 
[~^] = exp (-0. 33) ~ 0. 72 


1/2 
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Thus, approximately 28 percent of the incident radiation at 6. 3 pi will be absorbed 
by water vapor along a tangential path, whereas carbon dioxide, the major atmospheric 
constituent, is transparent at this wavelength. 


Attenuation by Scattering 

Attenuation of radiation at 6. 3 pi other than by water vapor absorption should be due 
only to scattering. Accordingly, Ben-Aryeh (ref. 4) has suggested that the resulting 
signal reduction be monitored with a second spectral channel at 8. 4 pi, where atmos- 
pheric absorption is expected to be negligible. 

Even under the most adverse condition, Rayleigh scattering can be shown to be 
negligible. The Rayleigh scattering cross section a of carbon dioxide at 0. 5890 pi is 
9. 18 x 10-28 cm 2 ( re f. 7). Assuming it varies as the inverse fourth power of the wave- 
length, it decreases to 7 x 10 “82 cm 2 a t 6. 3pt. Therefore, even in traversing a tangen- 
tial path, where the quantity of carbon dioxide is a maximum, i. e. , 9. 2 x 10^4 mol/ 
cm2-column [equation (13) below] , the fractional decrease of intensity by Rayleigh 
scattering is given by (ref. 8): 

_QO O A n 

I = oN = 7 x 10 x 9.2 x 10 = 6.45 x 10~' , 

s z, max 

a quantity too small to affect the proposed measurement. Loss by Mie scattering 
therefore appears to be the most likely mechanism of attenuation. It might be caused, 
for example, by dust clouds (ref. 8) or by frozen carbon dioxide particles, which have 
been suggested as an almost permanent feature of the Martian atmosphere (ref. 9). 
Specifically, therefore, a scattering measurement in the infrared would indicate the 
presence of particulate matter in the Martian atmosphere. 


Absorption by Carbon Dioxide 

From equation (9a), for carbon dioxide: 


N 

z, max 


2. 5 x 2 x 10 17 ( 336 O x 10 x 10 10 ) 
24 2 

9.2 x 10 mol/cm -column 


\ 1/2 


(13) 
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Since 1 cm at STP contains 2. 68 x 10 mol, the above amount of carbon dioxide cor- 
responds to 


9 - 2xl ° = 3.34 x 10 5 atm-cm. 

2.68 x 10 19 


(13a) 


Infrared absorption of carbon dioxide in terms of atm-cm is tabulated in the literature 
(ref. 7). 

Measurement of attenuation at 4. 3 ju will permit an estimate of the carbon dioxide 
profile in the Martian upper atmosphere. Extinction occurs when the total amount of 
carbon dioxide is slightly less than 2000 atm-cm (ref. 6). This quantity is exceeded by 
more than two orders of magnitude along the tangential path as deduced in equation (13a). 

2 

The extinction altitude (total number of mol/cm -column along path tangential to 
surface equal to 2000 atm-cm) is readily calculated from equation (9), i. e. : 

h o = H CO log e [“A ' 1 < 2 >' KH ) 1/2 ] 

= 10 log e |(^2 x 10^/2 x 2. 68 x 10^) ( 2n x 3360 x 10 x 10 

= 51 km . 



INSTRUMENTATION 

The experiment requires the spacecraft to travel on a trajectory behind the dark 
side of Mars and for the infrared sensors to point at the Sun. We assume a flyby dis- 
tance of 10, 000 km and a spacecraft velocity of 10 km/sec. 

Since, as discussed previously, a measurable indication of infrared absorption by 
water vapor cannot be expected beyond 20 km from the Martian surface, an altitude 
resolution of at least 1 km is desirable. At 10,000 km, this requires a field of view of 
20 by 20 sec. 

To show that adequate signal strength can be obtained with existing state-of-the- 
art instrumentation, the most adverse case of optical design will be considered. This 
consists of a 0. 1 by 0. 1 mm thermistor flake and a collimator with a field of view of 
20 by 20 sec, which simply consists of two 0. 1 by 0. 1 mm slits, 1 meter apart, with 
a suitable infrared filter. 
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Water Vapor Measurement 

Since the average distance between the Sun and Mars is 1. 52 A. U. , the value of 
the solar constant at Mars' distance is 0. 06 w/cm^. 

The spectral region between 5. 5 and 6. 5^i contains 0. 17 percent (ref. 10) of the 
total irradiance, i. e. , 1. 02 x 10 “4 w/cm^. Assuming a filter efficiency of 35 percent, 
the radiative power passing through the collimator is 3. 57 x 10 w. With a typical 
solid-backed bolometer infrared detector, for a sensitive area of 0. 1 by 0. 1 mm, the 
sensitivity is 2500 V/w. The resulting detector output will be 8. 95 jxV, whereas the 
Johnson noise of the bolometer resistors and the noise level of the preamplifier can be 
kept below 1 fj,V (ref. 11). The predicted 28 percent change of the incident radiation 
caused by absorption by water vapor will produce a signal of approximately 2. 5 fiV. 
Thus, a detection scheme of this type will be capable of establishing an upper limit of 
10~3 g/cm2-column of water vapor in the Martian atmosphere. This figure might well 
be lowered by a more efficient optical design (ref. 12) or by relaxing the altitude res- 
olution of the radiometer. 


Scattering Experiment 

The spectral region between 8 and 9 ju contains less than 0. 05 percent (ref. 10) of 
the total irradiance. Following the above procedure, a 28 percent change of incident 
radiation under identical conditions would not be detectable. On the other hand, the 
spectral region between 3. 4 and 3.6 u contains 0. 22 percent (ref. 10) of the total 
irradiance. In this case the unattenuated signal will be 12 MV and a change of less than 
20 percent will be readily detectable. 


Carbon Dioxide Measurement 

The spectral region between 4. 1 and 4. 5 fx contains 0. 2 percent (ref. 10) of the 
total irradiance. Under identical conditions of collimation, transmittance, and de- 
tector response, the unattenuated signal will be 10. 5 juV. Because of the high abundance 
of carbon dioxide in the Martian atmosphere , the carbon dioxide profile will be readily 
measurable. 

The basis of the instrumented package will be a radiometer design developed for 
the Mariner II mission (ref. 12) with an optical system considerably more sophisti- 
cated than the one outlined above. Except for certain modifications, such as the 
addition of a third channel and filters for different portions of the infrared spectrum, 
the equipment is already fully flight-qualified. 

A comparison of a Martian radiometer with that flown on Mariner II is given in 
Table I. 
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TABLE I 



Mariner H 

Mars 

Channel A: spectral range 

00 

CD 

5. 5 - 6. 5 fl (Water Vapor) 

Channel B: spectral range 

10.4 n 

3.4-3 .6 11 (Scattering) 

Channel C : spectral range 

- 

4. 1 - 4. 5 [i (Carbon Dioxide) 

Channel B: amplifier 

logarithmic 

all channels linear 

Field of View 

1. 2° x 1.2° 

20 x 20 sec 

Chopper Frequency 

20 Hz 

100 Hz 

Power at 400 Hz 

1. 9 w 

6 w 

Reference Optics 

essential 

redundant 

Reference Amplifier 

essential 

redundant 


Except for the chopper frequency, the above changes do not require further explan- 
ation. The figure for the chopper frequency is derived as follows: Assuming the most 
adverse case, the scan-speed will be equal to the spacecraft velocity, i. e. , 10 km/sec. 
(A reduction of the scan speed is, of course, possible by a suitable change in orientation 
of the velocity vector of the spacecraft relative to the Sun- Mars position vector. ) The 
time constant of the bolometer is 3 msec; for 99 percent response, at least five time 
constants are required, namely, 15 msec. For a 50 percent duty cycle, 30 readings 
per second can thus be obtained. The chopper frequency should be at least three times 
the readout frequency, i. e. , approximately 100 Hz. This can be obtained with a modest 
increase in chopper power from 1. 9 to approximately 4 w. 

Readings can thus be obtained every 0. 33 km. Regarding data acquisition, for an 
accuracy of 1 percent full-scale, seven binary bits are required; the bit rate is 210 per 
second per channel. Scanning from 100 km to the planetary surface requires a bit- 
storage capacity of 2100 per channel, i. e. , a total of 6300. This number is doubled if 
the scan can be performed on both sides of the Martian disc. 
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CONCLUSIONS 


The underlying theory of the proposed experiment relies on the best presently 
available data on the Martian atmosphere. No unreasonable assumptions have been 
made in the calculations. In essence , the instrumentation has been qualified on 
Mariner II (ref. 13) and the proposed modifications will not affect performance of 
reliability. Where necessary, the equipment will have to be updated in the light of 
advances made since 1962. 

Successful functioning of the equipment can be anticipated with confidence with the 
following tentative characteristics of the radiometer: 

Size: 6 by 6 by 6 inches 

Weight: 5 lb 

Power : 0. 5 w - 12 V dc , 6 w - 400 Hz , 26 V rms. 


National Aeronautics and Space Administration 
Electronics Research Center 
Cambridge, Massachusetts, April 1967 
185-47-26-01-25 
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